o 



X 



MAN/HEP/2010/5 
IPPP/10/36 
DCPT/10/72 



Extracting Higgs boson couplings using a jet veto 



B.E. Cox and J.R. Forshaw 
School of Physics & Astronomy, University of Manchester, Oxford Road. Manchester MIS 9PL. U.K. 

A.D. Pilkington 

Institute for Particle Physics Phenomenology, University of Durham, South Road. Durham DHl 3LE. U.K. and 
School of Physics & Astronomy, University of Manchester, Oxford Road. Manchester M13 9PL. U.K. 

- - - We show that the Higgs boson's effective couplings to gluons and to weak vector bosons can be 

' extracted simultaneously from an analysis of Higgs plus dijet events by studying the dependence of 

^ _ ^ , the observed cross-section upon a third-jet veto. 

(N : 

JIh '. 1. INTRODUCTION 

^ ■ 

I ■ If electroweak symmetry is broken by the Higgs mechanism, it is very likely that one or more Higgs bosons will be 
\ detected at the LHC. Measm^ements of the spin, CP and couplings of the observed Higgs states will be required to 
I— i. determine the nature of the Higgs sector. In this letter, we propose to extract the Higgs boson's effective-couplings to 
'■ weak vector bosons (Ay) and to gluons (Ag) by fitting the observed cross-section for pp — > Hjj + A as a function of a 
central jet veto scale, Qo- Specifically, we focus upon collisions in which a Higgs boson [H) is produced in between two 
O i' hard jets (j) that are far apart in pseudo-rapidity. Our interest is in studying the behaviour of these events subject to 
O , the further constraint that there should be no third jet with transverse momentum above Qo in the pseudo-rapidity 
1-^^' region between the two hard jets. 

' In the standard Higgs-plus- two-jet analyses, the central jet veto enhances the contribution of vector boson fusion 
Ps) (VBF) over gluon fusion (GF) by virtue of the difference in i-channel colour flow. In particular, Qq is usually chosen 
J> to be as low as possible whilst remaining robust against uncertainties due to the underlying event (UE). As Qo is 
, increased, the relative contribution oi gg ^ H is enhanced. Therefore, if the cross-section is measured as a function 
' of Qo, then a fit of the form a{Qo) = AgCr™((5o) + ^V^y^iQo) allows for the simultaneous extraction of the Higgs 
. couplings to both vector bosons and to gluons provided the Standard Model (SM) predictions (crg^((5o) and cry^(Qo)) 
are sufficiently well known. ^ 

To illustrate the point, we consider a 120 GeV Higgs boson decaying to r+r", although the analysis should also be 
' applicable to other decay channels and different Higgs masses. The parameters Ag and Ay measure deviations from 
the SM expectation and are sensitive to the following combination of Higgs boson decay widths: 



SMt3T}/tj , _ rggFT-.^ , . SU-r,-r,/TT . ^WW^T7 



Aga^^^ BRiH ^ rr) oc and Aya^^^^ BR{H ^ rr) oc . (1) 

^ ' I total I total 



The second relation assumes the ratio of WWH and ZZH couplings is the same as for the Standard Model, which is a 
generic expectation of custodially symmetric Higgs sectors. If new physics breaks this universality then Ay measures 
the effective coupling to weak vector bosons. By fitting to a linear combination of the SM cross-sections, we are 
assuming that any new physics generates deviations in the Higgs boson couplings but does not otherwise influence its 
production and decay.^ We note that the ratio Fgg/rww may be obtained directly from such a measurement. 



THE SIGNAL CROSS-SECTION 



We use the Sherpa 1.2 Monte Carlo (MC) [l| to generate sarnples of VBF and GF events in pp collisions at 
■^s = 14 TeV using the CTEQ6L parton distribution functions We generate H + n parton matrix elements 



^ The interference between the two is negligible. 

^ An example where this is not the case is the MSSM scenario where the Higgs boson is produced as a result of bottom quark fusion. 
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FIG. 1: The Vector Boson Fusion and Gluon Fusion contributions to the cross-section in the kinematic range described in the 
text for the process pp —i- H + jj + X with rr as a function of Qo , the central jet veto scale: Figure (a) shows the absolute 

cross-section in fb; Figure (b) shows the cross-sections normalized to unity at Qo = 50 GeV. The default samples are labelled 
as H + 3j (Qckkw = 15 GeV). Also shown are the effects of adding underlying event (UE) and changing the CKKW matching 
scale to Qckkw = 30 GeV. 



(n = 2,3) with CKKW matching between the matrix elements and the parton shower Q. We also invoke if- factors 
to ensure that the generated cross-sections, after matrix element plus parton shower matching, are equal to the 
corresponding NLO results of Campbell, Ellis & Zanderighi after implementing their 'weak boson fusion search cuts' 
We find that i^vBF = 1-18 and Kqf — 2.11. The latter re-scaling is comparable to the factor of 1.84 used by 
Andersen, Campbell and Hoche Q to match Sherpa to MCFM |6j] at y^s = 10 TeV. 

In order to simulate the experimental acceptance, we apply the same cuts as those used in a recent ATLAS study 
. The anti-fcx algorithm is used with radius parameter R — 0.4 to identify jets and jets that fall within a cone 
of radius AR = 0.2 around the true r direction are identified as coming from a r-decay.^ The two r candidates are 
required not to be back-to-back in azimuth, i.e. cosjAcj)) > —0.9, to facilitate the determination of the rr invariant 
mass using the collinear approximation. Following [7|, we apply further cuts on the highest transverse energy (Et) 
jets in the event excluding the r jets: 

Et,i > 40 GeV, Et^2 > 20 GeV, Mjj > 700 GeV, Arj > 4.4, r/i x r/2 < , (2) 

where r]i_2 are the pseudo-rapidities of the jets, A77 = \rii — 772] and Mjj is the dijet invariant mass. In addition, the 
missing transverse energy is required to be greater than 30 GeV. 

Figure [TJa) shows the Vector Boson Fusion and Gluon Fusion contributions to the cross-sections (with the X- factors 
applied and after cuts) as a function of the central jet veto scale Qo- FigurelTJb) shows the cross-sections as a function 
of Qo normalized to unity at Qo = 50 GeV in order to show the Qo dependence. Also shown are the effects of changing 
the CKKW matching scale from 15 GeV to 30 GeV and of adding the default underlying event in Sherpa. We place 
both of these uncertainties in context in the next section. 

In order to estimate the likely experimental event rates, the cross-sections shown in Figure [T] must be corrected to 
account for the experimental efficiencies ey and eg for selecting VBF and GF events respectively. These arise from 
the trigger efficiency, jet/r/lepton reconstruction and r/lepton identification. It is not possible to account for these 
effects in a study such as this, and we estimate ey = 0.036 for VBF using the numbers presented in Q for all r decay 
channels. The efficiency also includes a veto on ^-tagged leading jets to reduce the tt background. This is the only 
efficiency that is expected to vary between VBF and GF. However, this difference will be small relative to the overall 
normalization uncertainty on the GF cross-section and we therefore assume eg = ey = e. 

We extract Ag and Ay using a pseudo-experiment approach. The number of events for each process (VBF or GF) 
is determined using a Poisson distribution. The mean is given by eAiCTi(50 GeV)L, where L is the luminosity and ai 
is the cross-section (shown in Figure [ija)) at Qo — 50 GeV.^ Events are then selected at random from MC samples 



^ Experimentally, t's are identified using dedicated algorithms and jets that overlap with these are identified in the same way. 
^ At this value of Qo we expect 91 VBF events and 25 GF events after all cuts with 60 fb~^ of data. 




that survive the experimental cuts detailed above and a further cut on Qo allows the Qo distribution to be generated. 
Each point in the Qo distribution can then be smeared according to the expected systematic uncertainties, discussed 
below. The GF and VBF distributions are then combined to create the pseudo-data set. 

To extract the couplings, we fit the pseudo-data set to a function of the form f{Qo) = Ag/g(Qo) + Av/v(Qo), where 
fiiQo) — ^Lai {Qo) are the SM theoretical predictions (i.e. the solid curves in Figure [ija)). The whole procedure is 
repeated 10"* times for each value of Ag and Ay and the variation in the fitted values of Ag and Ay determines the 
accuracy to which we expect to measure Ag and Ay. 

Figure [DJa) shows the expected uncertainty in the extracted value of Ag as a function of both Ag and Ay for 
L — fb~^. The scale on the right of the figure refers to the fractional uncertainty, JAg^y/Ag^y. Figure[2Kb) shows 
the corresponding uncertainty in Ay. At this stage, we have ignored all theoretical uncertainties and experimental 
systematic errors. That Ay should be determined with greater accuracy over much of the parameter space is a result 
of the fact that the VBF cross-section is typically larger than the GF cross-section after cuts. 



3. EFFECT OF SYSTEMATIC UNCERTAINTIES 



The dependence of the cross-section on the central jet veto, ct(Qo)i can be expressed as 

'^(Qo) =^Tjj(l--Pvcto(Qo)) (3) 

where CTjj is the Hjj cross-section with no jet veto and Pvcto(Qo) is the probability of finding a third jet above Qo 
in the inter-jet region. Uncertainties in ajj affect the overall normalization whilst uncertainties in 1 — Pvcto(Qo) also 
affect the shape of the Qo distribution. We can account for the theoretical uncertainty in the shape by estimating the 
uncertainties in the predictions for Pvoto at Qq = 20 GeV and at Qq = 50 GeV, which we treat as uncorrelated. In 
practice this means shifting the data points at these two extremal values of Qo in each pseudo-data set according to 
the theoretical uncertainties. For intermediate values of Qo we shift the pseudo-data points by interpolating linearly 
in Qq between the extremal shifts. 

Uncertainties in the VBF cross-section as a function of a third jet veto have been studied explicitly in the literature 
0i estimate of the dominant NLO corrections to Hjjj were studied, and it was concluded that PvctolQo) 
is known to an absolute accuracy of better than 1% for all relevant values of Qq. In [1^], the NLO Hjj cross-section 
was matched to the herwig and pythia parton showers in order to estimate -Pvcto(Qo)- The uncertainty was found 
to be larger (±3%) and comparable to the result of [oj for the LO Hjjj calculation. We assume here that the 
yet-to-be-calculated NNLO Hjj calculation with parton shower matching will confirm the NLO Hjjj result of @. 

The overall normalization of the VBF cross-section is the other main theoretical uncertainty. The partial NNLO 
calculation presented in concludes that the overall uncertainty due to unknown higher-order QCD corrections 
is around 2% (at y^s = 7 TeV and without any VBF cuts). A complete NLO calculation, including electroweak 
corrections, is presented in [l^l where it is concluded that the total Hjj cross-section with VBF cuts should be known 
to better than 2% for Higgs masses in the range 100-200 GeV. On top of all these corrections, we must account 
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for an overall uncertainty of some 3% arising from uncertainties in the parton distribution functions [13l Il4| . The 
net effect is that perturbative uncertainties on the VBF cross-section are probably no more than ±4% in the overall 
normalization and with no significant error on the Qq dependence. We implement this normalization uncertainty by 
varying all pseudo-data points by ±4%, according to a fiat distribution in Qo- 

The theoretical uncertainties on the GF cross-section are considerably larger and the accuracy to which the Qo 
dependence of the cross-section can be computed is not well known. The uncertainty in the NLO Hjj calculation in 
the absence of a veto (o-jj) was estimated in [J] to be around ±20% and we assume this value. We assign an additional, 
uncorrelated, uncertainty of ±20% to 1 — Pveto(Qo) at Qo = 20 and 50 GeV. The study in suggests that the theory 
may be rather more uncertain at the present time [H, \ % [l^ . However, an uncertainty of the order that we assume 
here should be attainable in the near future (see also [l6| for a study of the Qo-dependence induced by wide-angle 
soft-gluon emission). The uncertainty due to CKKW matching illustrated in Figured] is clearly small on the scale of 
these other uncertainties. 

Perturbative uncertainties are not the entire story. Chief amongst other theory uncertainties is the lack of knowledge 
of the underlying event (UE). In Figure [U we illustrate the potential impact of the UE. The blue dot-dashed curve 
shows our prediction after including Sherpa's simulation of the UE and the blue solid curve is that in the absence of 
any UE. We assume an uncorrelated uncertainty of ±1% in 1 — Pvoto(Qo) at both Qo = 20 GeV and 50 GeV in the 
VBF cross-section. In the case of GF, the effect of the UE is increased and we assign an uncorrelated ±3% uncertainty 
at Qo = 20 GeV and 50 GeV. Early LHC data will allow tuning of UE models which should reduce the degree of 
uncertainty. 

In addition to theoretical uncertainties we need also to estimate the impact of experimental systematics. The 
dominant systematic is due to the uncertainty in the jet energy scale (JES). The ATLAS collaboration estimated the 
total systematic uncertainty on the VBF cross-section to be ±20% at Qo = 20 GeV Q, of which 16% is due to the 
JES and 10% to other sources. We can satisfactorily reproduce the uncertainty arising from the JES by shifting the 
generator level jet energies by the ATLAS defaults (i.e. 7% for jryj < 3.2, 15% otherwise). Importantly, we observe no 
significant dependence of the uncertainty from the JES on Qo and therefore assign an overall normalization uncertainty 
to the VBF cross-section of ±20%. The uncertainty from the JES is larger for GF events, due to the steeper leading jet 
i?T and Afjj distributions, and it has a mild Qo dependence. We therefore estimate an overall normalization uncertainty 
of ±30% in conjunction with an additional uncorrelated uncertainty of ±3% in 1 — Pvcto(Qo) at Qo = 20 GcV and 
50 GeV. We also assume that the effects from the pile-up of multiple proton-proton interactions will be well understood 
and do not assign an additional uncertainty from this source. 

To account for the influence of background processes, we assume that statistical fluctuations dominate.*^ This should 
be the case if the background can be extracted by fitting the distribution in data. We assume equal numbers 
of signal and background events at Qo = 50 GeV after cuts, which is broadly in line with the conclusions in Q- We 
also assume that the background varies with Qo such as to lie halfway between the upper and lower curves in Figure 
[TJb). Our treatment of the background is clearly very approximate but is sufficient for the purpose of generating 
additional statistical fluctuations across the Qo distribution. We note that the Qo dependence of the background will 
be determined from data. 

We are now in a position to quantify the effect of these various systematic uncertainties. TableUillustrates the effect 
of the individual experimental and theoretical uncertainties on the extraction of the effective couplings and Figure [3] 
shows the expected uncertainty after the inclusion of all systematic uncertainties. In the table, the numbers in the 
first row correspond to ignoring all experimental and theoretical uncertainties except those due to statistics and those 
in the final row represent the uncertainties after accounting for the various uncertainties discussed above. The middle 
four rows show the effect of including only the uncertainties on (i) statistical fluctuations in the background; (ii) the 
VBF cross-section; (iii) the GF cross-section; (iv) experimental systematics. We also show results for a statistically 
comparable 'BSM' point (the couplings are chosen so that the cross-section is approximately equal to that in the SM 
at Qo = 50 GeV). 

4. CONCLUSIONS 

The method we propose permits the measurement of the ratio of the effective couplings of the Higgs boson to gluons 
and to weak vector bosons within a single analysis. If the branching ratio to t leptons is known, then the method 
can be used to extract the individual couplings. The couplings can be determined to an accuracy that is comparable 



^ The backgrounds are discussed in detail in Zjj is the most important background but tt production, W^-|-jcts production and QCD 
multi-jet production are also relevant. 
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FIG. 3: The uncertainty in the extracted values of (a) Ag and (b) Av given the default ATLAS cuts (see text) and 60 fb"'^ of 
data. The effects of systematic errors are estimated as discussed in the text. The scale on the right of each figure refers to the 
fractional uncertainty, i.e. (5Ag.v/Ag,v- 



Error 


SM (Ag 

CTAg/Ag 


,v = l) 

(TAy / Av 


BSM (Ag = 

cr Ag/Ag 


= 4, Av = 1/4) 

(TAv / Av 


Stat, only 


0.51 [0.23] 


0.16 [0.07] 


0.19 [0.08] 


0.72 [0.33] 


Backgd. 
VBF 
GF 
Expt. 


0.56 [0.25] 
0.52 [0.25] 
0.65 [0.45] 
0.62 [0.39] 


0.18 [0.08] 
0.17 [0.08] 
0.19 [0.11] 
0.26 [0.21] 


0.20 [0.09] 
0.19 [0.08] 
0.43 [0.40] 
0.35 [0.31] 


0.79 [0.35] 
0.75 [0.33] 
1.56 [1.40] 
0.89 [0.52] 


AU 


0.77 [0.57] 


0.28 [0.23] 


0.53 [0.50] 


1.66 [1.49] 



TABLE L Fractional error in extraction of Ag and Av for the SM and a statistically comparable BSM parameter point. Numbers 
correspond to 60 fb^^ of data or, in square brackets, to 300 fb^^ of data. 



to other methods [13, 13 • A model- independent measurement of the relative size of the VBF and GF contributions 
should be of considerable value in subsequent analyses to study the CP nature of the HWW pj|] and Htt [20, HH 
couplings via the azimuthal angle dependence of the tag-jets. 

Our results demonstrate the need for an improved theoretical understanding of the veto dependence of the gluon 
fusion contribution. The primary experimental uncertainty comes from the knovi^ledge of the jet energy scale. 

We have not made any attempt to optimize the experimental cuts for this analysis, preferring instead to follow 
the study presented in [7|, which was optimized for VBF production. It should therefore be possible to improve on 
our estimates of the coupling to gluons. Our method could be further generalized to include other Hjj production 
channels, e.g. bb H in the case of MSSM scenarios in which the Higgs has an enhanced coupling to 6-quarks. 

Using early LHC data, important precursors to the analysis presented here would be to perform analyses of the 
jet-veto dependence in dijet [22|, W/Z + j and W/Z + jj events (23| . 
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